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Magnetic properties of a cr— Fei6Cri4 alloy calculated with the charge and spin self- consistent 
Korringa-Kohn-Rostoker (KKR) and combined with coherent potential approximation (KKR-CPA) 
methods are reported. Non-magnetic state as well as various magnetic orderings were considered, 
i.e. ferromagnetic (FM) and more complex anti-parallel (called APM) arrangements for selected 
sublattices, as follows from the symmetry analysis. It has been shown that the Stoner criterion 
applied to non-magnetic density of states at the Fermi energy, Ef is satisfied for Fe atoms situated 
on all five lattice sites, while it is not fulfilled for all Cr atoms. In FM and APM states, the values 
of magnetic moments on Fe atoms occupying various sites are dispersed between and 2.5 /j,b, and 
they are proportional to the number of Fe atoms in the nearest-neighbor shell. Magnetic moments 
of Cr atoms havin much smaller values were found to be coupled antiparallel to those of Fe atoms. 
The average value of the magnetic moment per atom was found to be < ^ >= 0.55/^s that is by 
a factor of 4 larger than the experimental value found for a a— Feo.53sCro.462 sample. Conversely, 
admitting an anti- parallel ordering (APM model) on atoms situated on C and D sites, according 
to the group theory and symmetry analysis results, yielded a substantial reduction of < /i > to 0.20 
\iB- Further diminution of < fi > to 0.15 /is, which is very close to the experimental value of 0.14 
fiB, has been achieved with the KKR-CPA calculations by considering a chemical disorder on sites 
B, C and D. 



I. INTRODUCTION 

Sigma (cr) phase has a tetragonal unit cell (type 
P42/mnm) hosting 30 atoms distributed over five crys- 
tallographically non-equivalent sites. Due to high coor- 
dination numbers (12-15) the cr-phase is a member of 
the Frank-Kasper family of phases. The importance of 
this class of structures is reinforced by the fact that they 
exhibit topological properties similar to simple metallic 
glasses due to an icosahedral local arrangement. Conse- 
quently, they can also be regarded as very good approxi- 
mants for dodecagonal quasicrystalsi. The cr— phase has 
been revealed to occur only in such alloy systems in which 
at least one element is a transition metal. 

Concerning the existence of the sigma-phase in bi- 
nary alloys, about 50 examples were found so far— 
Its physical properties are characteristic of a given al- 
loy, a high hardness and brittleness being a common 
feature. In the family of the binary alloy sigma-phases 
only two viz. Fe-Cr and Fe-V have well-evidenced mag- 
netic properties^—. In particular, the magnetic prop- 
erties of a— phase strongly depend on composition. In 
the Fei-^Cra; system, where the cr— phase exist in a nar- 
row range of ~ 0.45 > x >~ 0.50, the Curie temper- 
ature varies between K (for x w 0.50) and 38 K (for 
x « 0.45) 6 . In the Fe-V system, where the range of the 
(T-phase occurrence is about six-fold wider, the effect of 
the composition is even stronger. 

Another common characteristics is a lack of stoichiom- 
etry, a feature that in combination with the complex crys- 
tallographic structure, makes the cr-phase a very chal- 
lenging object both for theoretical investigations and in- 
terpretation of experimental data. On the other hand, 



the lack of stoichiometry may also have some advantage 
viz. the cr-phase can be formed in a certain range of 
composition, giving thereby a chance for studying the 
influence of content on its properties, hence a possibility 
of tailoring these properties. For example, a magnetism 
of the cr-phase in Fei-^Va, alloys is strongly dependent on 
the composition, and, in particular, the Curie tempera- 
ture may vary between K (for x « 0.65) and 320 K (for 
x 0.34) 8 . The complex structure and the chemical dis- 
order (all five sites are occupied by various atoms) result 
in a huge number of possible atomic configurations which 
makes an interpretation of measurements performed by 
using microscopic methods e.g. Nuclear Magnetic Res- 
onance or Mossbauer Spectroscopy to be a difficult and 
often non-unique task. In these circumstances, theoreti- 
cal calculations, as these presented elsewhere*^, can be 
very helpful as they not only provide a deeper insight into 
the underlying physics but they also offer a very useful 
support for a correct interpretation of experimental data. 

This paper reports theoretical calculations on the mag- 
netic properties of the cr— FeCr system which is regarded 
as the archetype of all cr— phases. Its magnetic proper- 
ties were recently intensively investigated^, and, in par- 
ticular, the average magnetic moment per Fe atom was 
determined from magnetization measurements assuming 
only Fe atoms having parallel moments contribute. How- 
ever, the knowledge is not complete yet, as it is unknown 
whether or not both types of atoms (Fe and Cr) are mag- 
netic, and if so, what are the values of their magnetic mo- 
ments on particular sublattices and what is their mutual 
orientation. These questions cannot be answered exper- 
imentally because one cannot grow big enough single- 
crystals of the sigma for a decisive neutron diffraction 



2 



experiment. Hence, the situation prompts and justifies 
further relevant theoretical calculations. The ones re- 
ported in this paper are in the line with this aim. 




FIG. 1: (Online color) Unit cell of the a-phase. Atoms be- 
longing to various crystallographic positions are indicated by 
different fillings. NN atoms of each of five lattice sites are 
indicated also. 



II. COMPUTATIONAL DETAILS 

As aforementioned, the cr-phase has complex close- 
packed tetragonal structure with thirty atoms in the 
unit cell (Fig. [T]). Atoms are distributed over five non- 
equivalent sites in the unit cell, called A, B, C, D and E, 
the population of them is shown in Table HI Each posi- 
tion can be characterized by the total number of nearest 
neighbors, NN, their distribution over five sublattices, 
the distances to NN atoms as well as the occupancy by 
Fe or Cr atoms. Some properties can be derived directly 
from the space group information. Each atom in the er- 
phase structure has a high coordination number (from 
12 to 15 atoms), belonging to various sites (Table HJ. 
Their average inter-atomic distances substantially differ 
and are in the range from ~ 2.27 A(the smallest E-E 
distance) to ~ 2.92 A(the largest B-E distance). In- 
terestingly, the NN spatial distribution is not far from 
spherical for each site (Fig. [TJ. 

Electronic structure calculations for cr— FeCr com- 
pounds have been performed using the charge and spin 
self-consistent Korringa-Kohn-Rostoker method. The 
crystal potential of the muffin-tin form has been con- 



TABLE I: Atomic crystallographic positions and numbers of 
the nearest neighbor atoms, NN, for the five lattice sites of 
the cr-phase. 



Site 


Crystallographic positions 


NN 






A 


B 


C 


D 


E 


Total 


A 


2i (0, 0, ) 




4 




4 


4 


12 


B 


4f (0.4, 0.4, ) 


2 


1 


2 


4 


6 


15 


C 


8i (0.74, 0.66, ) 




1 


5 


4 


4 


14 


D 


8i (0.464, 0.131, ) 


1 


2 


4 


1 


4 


12 


E 


8j (0.183, 0.183, 0.252) 


1 


3 


4 


4 


2 


14 



structed within the local density approximation (LDA) 
framework, applying the von Barth-Hcdin formula^ for 
the exchange-correlation part. With the self-consistent 
crystal potentials converged up to 0.1 mRy and charges 
up to 10~ 3 e, total, site-decomposed and ^-decomposed 
(s, p and d) density of states (DOS) have been com- 
puted on a 601 energy point mesh for ordered models, 
employing the tetrahedron k-space integration technique 
(dividing the irreducible part of the Brillouin zone into 
120 small tetrahedrons). 

In the first stage, the electronic structure calculations 
have been done for ordered approximants of the cr— FeCr 
alloys, and for that purpose we lowered the symmetry of 
the unit cell (space group P42/mnm) to a simple tetrag- 
onal one. In practice, the tetragonal unit cell and atomic 
positions were unchanged but variable occupancy made 
all 30 atomic positions crystallographically nonequiva- 
lent. In such a specified unit cell each of the crystallo- 
graphic positions was occupied exclusively either by Fe 
or Cr atom. However, in our numerical runs we were con- 
strained by the experimentally determined Fe/Cr concen- 
trations on each of the five lattice sites^ and the com- 
position of the cr— FeieCri4, being close to the measured 
stoichiometry (Feo.53sCro.462), has been considered. 

It should be noticed that the composition of the al- 
loy is given in two equivalent notations. When the unit 
cell of the cr-phase is referred to calculations, we use the 
Feso-ajCrj, formula, since x must be the whole number 
here. Conversely, the Fei-xCra, formula is related to 
those cases, when x (the Cr concentration) in the sample 
was obtained experimentally. 

Then, the effects of chemical disorder on the electronic 
structure have been studied employing the KKR method 
with the coherent potential approximation (CPA), which 
belongs to well-established techniques used in a descrip- 
tion of intermetallic alloys. In this case, the space group 
P42/mnm, defining unit cell of cr— phase, was allowed 
and five nonequivalent sites were occupied as follows: A 
(2i site: 100% Fe), B (4/ site: 25%Fe and 75%Cr), C 
(8i site: 25%Fe and 75%Cr), D (8i site: 100%Fe) and 
E (8j site: 37.5%Fe and 62.5%Cr). The occupancies of 
the sublattices intentionally differed slightly from those 
used for the ordered approximants (although it is easy to 
apply precise concentrations in the KKR-CPA method) 
due to the fact that we intended to make possible a coher- 
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TABLE II: Fe site-occupation parameters of the cr-FeCr alloy, 
experimental and assumed for calculations. N t stays for the 
percentage of the total Fe atoms on the site (referred to the 30 
atoms unit cell), whereas N describes the number of Fe atoms 
occupying the site. Ni and N2 represent the values used in the 
calculations for ordered (KKR) and disordered (KKR-CPA) 
models, respectively. The corresponding relative occupancies 
(percentage) are given in parenthesis. 



Site 


N t 


N 


Ni 


iV 2 


A 


11.3 


1.826 (91.3) 


2 (100.0) 


2 (100.0) 


B 


6.4 


1.040 (26.0) 


1 (25.0) 


1 (25.0) 


C 


20.5 


3.304 (41.3) 


3 (37.5) 


2 (25.0) 


D 


44.9 


7.208 (90.1) 


7 (87.5) 


8 (100.0) 


E 


17.0 


2.744 (34.3) 


3 (37.5) 


3 (37.5) 



ent comparison between ordered and disordered models. 
Thus, in the partly disordered model, Fe and Cr atoms 
were randomly distributed on three sites, while two other 
sites were constrained to be fully ordered. Consequently, 
this model slightly differs from the experimentally deter- 
mined distribution of Fe and Cr atoms in the unit cell 
of the a— FeCr phase. In real samples a small fraction of 
chromium atoms was also detected at A and D sites (see, 
Table HH) . However, we believe that the partly disordered 
model, forced by highly time-consuming KKR-CPA cal- 
culations, appears as a good approximation of a fully 
disordered cr— FeCr alloy. 

In KKR-CPA calculations, the spin-dependent charge 
and CPA cycles have been performed self-consistently 
on the complex energy plane using an elliptic contour 
(divided into 12 sections with 4 Gaussian quadrature 
points). The KKR-CPA Green function was computed 
on 75 special k point mesh in the irreducible part of the 
Brillouin zone. The Fermi energy, Ep, for all investigated 
alloys has been determined from the generalized Lloyd 
formula?^, which permitted us to find Ep precisely from 
the full derivative of the CPA Green function (without 
integrating over occupied states). Like for KKR runs on 
ordered models of the cr— FeCr alloy, for finally converged 
charges and crystal potentials in the aforementioned dis- 
ordered model, the KKR-CPA calculations resulted in 
total, site-decomposed and I— decomposed (with l max —2) 
DOSs computed on a 251 energy point mesh. More de- 
tails on the KKR-CPA method can be found in Refs . 14 ' 15 . 

In all our calculations, experimental crystallographic 
data (lattice constants and atomic positions) were 
employed^. 

In order to study the magnetism of the cr— FeCr alloy 
spin-polarized KKR as well as KKR-CPA calculations 
have been performed. A collinear ordering of magnetic 
moments (on Fe and Cr atoms), as allowed by the em- 
ployed computational codes, were assumed in a ferromag- 
netic (FM) state both in the ordered approximants and 
in the disordered alloys. However, the FM model of the 
cr— FeCr alloy yielded a serious discrepancy between cal- 
culated and experimental values of a total magnetization. 
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FIG. 2: (Online color) Average magnetic moments, < fi >, 
for five crystallographic sites versus average < Id ■ nd(Ep) >- 
value. Full symbols stand for Cr atoms, while empty ones for 
Fe atoms. 

Consequently, a more refined search for achieving a better 
arrangement of Fe/Cr magnetic moments, with the use of 
the symmetry group analysis (see Sec. Ill C), have been 
carried out. A relatively complex model of the magnetic 
structure, based on anti-parallel alignments of the mag- 
netic moments within selected sublattices (called APM) , 
was proposed and applied. Finally, spin-polarized KKR- 
CPA calculations were also performed for such magnetic 
unit cell, assuming a collinear ordering of Fe and Cr mag- 
netic moments, only. 

In all figures presented throughout the paper, Ep is 
located at zero, and spin-polarized DOS curves are given 
in Ry _1 per spin direction. 

III. RESULTS AND DISCUSSION 
A. Nonmagnetic state 

We start the analysis of electronic and magnetic prop- 
erties of the a— FeCr from the spin non-polarized cal- 
culations in order to verify the origin of magnetism in 
view of the Stoner model. Table IIIII presents averaged 
computed values of d-DOS at Ep, rid(Ep), and corre- 
sponding Stoner parameter, Id, on all five nonequiva- 
lent sites occupied by Fe and/or Cr atoms. We see that 
the Stoner criterion (Id ■ rid(Ep) > 1) is satisfied for all 
iron atoms (I^ e ~ 0.036Ry), being the highest one in 
the case of Fe at B site and the smallest one for Fe at 
A site. Conversely, the same condition is not fulfilled 
for chromium atoms (l9 r ~ 0.034i?y) essentially due to 
much lower DOS at Ep, which would suggest that the 
magnetic moments appearing on chromium atoms in the 
cr— FeCr alloy are in principle induced due to the presence 
of rather large values of Fe magnetic moments within the 
NN shell. Figure shows a dependence of the averaged 
magnetic moments calculated for the five sublattices ver- 
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sus the Stoner product that well supports the aforemen- 
tioned experimental results, suggesting itinerant charac- 
ter of magnetism in this system 7 -. 

B. Ferromagnetic state 

The aforementioned predictions from the spin non- 
polarized DOS and the Stoner product analysis have been 
followed by spin-polarized KKR calculations in order to 
determine local and average magnetic moments on Fe 
and Cr atoms. The magnetic structure of the a— FeCr 
has not been determined yet. Hence, we first discuss the 
ferromagnetic ordering, being the expected upper limit 
of total magnetization. The values of the magnetic mo- 
ments of each of 30 atoms in the unit cell obtained for the 
ordered approximants are presented in Fig. [3] One should 
notice here a relatively large dispersion of the magnetic 
moment of Fe atoms, \iFe-, from almost to nearly 2.5/is. 
The values of fj,p e for each of the five sublattices are, ad- 
ditionally, nicely correlated with the number of Fe atoms 
in the first coordination shell, NNp e . The linear corre- 
lation was described using the equation 

Mi> e = oi + biNN Fe (1) 

where index I =A,B,C,D,E denotes the particular sub- 
lattice. 

The average value of /i/ Fc , < /i/ Fe >, for each sublat- 
tice depends both on aj and bj parameters obtained from 
Eq. [1] fitted to the data. Using probabilities of finding 
particular AWpe-values, Pj(NNp e ), determined from a 
binomial distribution for each sublattice, < /ij Fe > can 
be calculated with the simple formula 

< W,« >= Hi Fa Pi{NN Fe ) (2) 

As is clear from Fig. < /Lt/ Fe > almost linearly in- 
creases with the average distances to their nearest atomic 
neighbors, < d >, for all sublattices. The magnetic mo- 
ments on Cr atoms are much smaller and mostly oppo- 
sitely aligned to those of Fe atoms. Their values range 
from -1.1 to 0.6 [is, and their correlations with the num- 
ber of NNpe are much weaker than those in the case of 
Fe. For this reason, the average values of hi Ct for vari- 
ous sublattices were obtained as arithmetic means of all 
KKR calculated values. 

Since, according to the present calculations, both Fe 
and Cr atoms possess the magnetic moment, it seems 
more reasonable to express the calculated magnetization 
of the a - FeCr phase, namely the total magnetic moment 
of the unit cell, per one atom in the unit cell (i.e. divided 
by 30) rather than per Fe atom (i.e. divided by 16 in this 
case). 

The average value of the total magnetic moment of the 
unit cell can be determined in two ways (Table [III]). 

1. As the average value of just all total magnetic mo- 
ments obtained in each calculation (e.g. 26 different ap- 
proximants of the cr— phase, involving 780 atoms). In 



this case, we are able to take into account the spin polar- 
ization of electrons located outside the muffin-tin spheres, 
giving non- vanishing contribution to total magnetization. 
This procedure gives 0.55/^s for the magnetic moments 
inside the muffin-tin spheres and only — 0.02/is for the 
interstitial area (0.53/xb in total). 

2. As the sum of the weighted five < /// >-values, 
where the weights are the numbers of atoms belong- 
ing to each sublattice (Eq. [2]). This method takes into 
account all possible configurations within the NN-shell, 
even those that have not been considered in the applied 
model. Interestingly, the obtained value of 0.59/zb is 
quite close to the previous one. 

An illustrative example of a spatial distribution of the 
magnetic moments of individual atoms in the unit cell is 
shown in Fig. [4^,. In order to present the calculated mag- 
netic structure more clearly, the z-axis projected unit cell 
was doubled along x and y directions. It is worth noticing 
that large magnetic moments (polarized parallel) assem- 
ble to some extend. Other analyzed configurations of 
atoms (not shown) give qualitatively similar picture. 

Comparable results were obtained in the work of Pavlu 
at al^. The authors performed spin-polarized calcu- 
lations using the VASP program, considering only the 
fully occupied sublattices by one type of atoms. For 
a cr— Feo.533Cro.467 sample the total magnetic moments 
(per atom) found by these authors were as large as 0.52, 
0.59 and 0.73 \ib for three possible atomic configura- 
tions. Applying similar simplifications for atoms distri- 
bution in the unit cell, Kabliman el al^- obtained /i=0.75 
[i B (o— Fe .6oCr .4o) and 0.25^ B (cr-Feo.47Cro.53) using 
FLAPW method. 

Due to the fact that these author a 16 ' 17 reported the 
average values of the magnetic moments on each sub- 
lattice and also for all considered atomic configurations, 
their data can be directly compared with our results. As 
can be seen in Fig. [51 in most cases the KKR values of 
the magnetic moments are slightly larger than the cor- 
responding ones reported elsewher o 16 ' 17 . The observed 
difference probably results from a difference in lattice 
constants taken for calculations in the present work (ex- 
perimental data) , being slightly different than those used 
in Refs. [l6| and[l3 (obtained from a relaxation of the unit 
cell) . Indeed, applying the a — b and c values of the lat- 
tice constants from Ref. [HI yielded smaller magnetic mo- 
ments. So, the overall agreement with the previous calcu- 
lations is satisfying, except for the sublattice E, where the 
results reported in Ref. [16| remain practically constant 
( 1.5 fj,B for Fe), while the corresponding values obtained 
in the present paper are clearly N Np e — dependent. Ac- 
tually, we observe that the magnetic moments on all five 
sublattices vary with NNp e in more or less similar way. 
However, we are not aware of any particular reason for 
markedly different behavior of Fe and Cr magnetic mo- 
ments on the E site with respect to the behavior found 
for other sites, as suggested in Refs [l6| and [l?]. 

A good compatibility of the theoretical results ob- 
tained for the fully ordered systems of the cr— FeCr 
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TABLE III: The average Stoner product, Si =< /d - ji^-Ef) >, (i=Fe,Cr) in the er-FeCr alloy as calculated for the non-magnetic 
(NM) state. Magnetic moments in /is of Fe and Cr atoms as well as their average values, < fj, >, are indicated for various 
magnetic models and particular lattice sites. jLt_p ei and /iFe 2 m the FM model stand for the values determined in two ways 
described in the text. The lack of values for Cr atoms on certain sites is a consequence of the assumed distribution as presented 
in Table [U 



Site 


NM model 


FM model 


APM model 


APM-CPA model 




SFe 


Scr 






P-Cr 


< n > 


fJ-Fe 




< /x > 


1-lFe 


HOr 


< n > 


A 


0.98 




0.77 


0.76 




0.76 


0.39 




0.39 


0.32 




0.32 


B 


1.26 


0.82 


2.16 


2.13 


-0.27 


0.33 


2.07 


-0.11 


0.43 


2.08 


-0.10 


0.44 


C 


1.33 


0.73 


2.08 


1.91 


-0.13 


0.63 


0.07 


-0.06 


-0.01 


-2.06/+2.11 


-0.11/+0.01 


-0.07 


D 


0.98 


0.68 


0.91 


0.90 


-0.04 


0.78 


0.24 


-0.01 


0.20 


-0.44/+0.65 




0.21 


E 


1.20 


0.67 


1.49 


1.38 


-0.39 


0.28 


0.99 


-0.20 


0.25 


0.90 


-0.25 


0.18 


av. 






1.30 


1.24 


-0.25 


0.55 


0.49 


-0.12 


0.20 


0.45 


-0.16 


0.15 



phas o 16 ' 17 assumed throughout the whole range of con- 
centrations on the one hand, and those obtained for the 
ordered approximants of disordered system (correspond- 
ing to realistic configurations for the a— Feo.533Cro.467 
alloy) on the other, allow us to deduce the influence of 
the NNp e on the hyperfine parameters to be weakly (if 
at all) dependent on the alloy stoichiometry. Therefore, 
the determined dependences of ni(NNp e ) for each sub- 
lattice can be assumed to be valid in the whole range of 
the cr-FeCr existence. 

The calculated average value of the magnetic moment 
(0.59 hb per atom per unit cell) should be referred to 
the experimental value obtained for the a— Feo.53sCro.462j 
viz. [i = 0.14/i^, which is almost four times smaller 
than the theoretical result. Such a large and unusual dis- 
crepancy between calculated and measured values of the 
magnetization on one hand, as well as a non-saturated 
behavior of M(i5f)-dependence£ on the other, may indi- 
cate the magnetic structure of the tr-phase is not a simple 



ferromagnetic one as was assumed in the aforementioned 
KKR calculations. 



C. Symmetry analysis and models of magnetic 
ordering 

The large discrepancy between the theoretical (assum- 
ing FM model) and experimental results of the total 
magnetization in the cr-FeCr prompted us to search for 
a possible different magnetic ordering in such complex 
phase using the group theory and symmetry analysis as 
described below. 

Let us consider an arbitrary crystal property, X, that 
is defined locally as u, for all atomic positions in a crystal 
with n atoms per unit cell. The local single-site property 
can be a scalar, a vector (polar or axial) or even a tensor. 
According to the symmetry analysis method, the global 
quantity X can be written in the most conveniently se- 
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FIG. 4: (Online color) Magnetic moments of Fe and Cr atoms 
in the unit cell, as seen along z axis, for (a) FM structure and 
(b) APM (i.e. anti-parallel ordering on D and C sites) struc- 
ture. Full and empty circles represent magnetic moments par- 
allel (/%,) and antiparall (/i a ) to the z— direction, respectively. 
The size of red (Fe) and green (Cr) circles is proportional to 
the /x- value. 



lected basis, what usually denotes a symmetry-adapted 
basis, consisting of the basis vectors of irreducible 

representations (IR) rff} of the symmetry space group, 
with as linear combination coefficients^, 



X = 



u(f{) 
u{f 2 ) 

U(fi +t Pl q lS ) 
u{f* 2 +tp,q,s) 

^yfl tp,q,s ) 



a (rl) 



: pa + qb + sc 



*^ ) (r;)e ifelF 

(3) 
(4) 



where i, n and A correspond to the vectors k, to the 
irreducible representations of the group and to the IR 
dimensions, respectively. 

In the case of any phase transition, the form of the ba- 
sis vectors and the information on the representations rel- 
evant to this transition (so-called active representations) 
are directly given by the theory of groups and the repre- 
sentations. One is able to calculate these quantities us- 
ing a dedicated software^. Then the linear combination 
coefficients can be determined by imposing proper condi- 
tions on the solution (specific for a given physical quan- 
tity) , and the symmetry-adapted ordering modes can be 
found. The representation r„ and the coefficients C^ ,v 



TABLE IV: Active representations of all types of modes for 
P42/mnm space group and k=0. The meaning of the symbols 
is described in the text. 



site 


mode 






l 


epresentation 








type 


n 






Ti 


TS 


TV, 


T7 


TS 


T>) 


T 10 




Pi 


l 


_ 


_ 








1 




_ 


_ 


2a 


P 2 


_ 


_ 


_ 


1 


_ 


1 


_ 


_ 


_ 


2 




Pa 




_ 


1 


_ 


1 


_ 


_ 


_ 


2 


_ 




Pi 




_ 


_ 


_ 


_ 


_ 


1 


_ 


_ 


1 


4/ 


P 2 


I 


- 


1 


1 


1 


1 


1 




1 


2 




Pa 




1 


1 


1 


1 


1 




1 


2 


1 




Pi 


l 




1 




1 




1 






2 


8i 


P 2 


2 


1 


2 


1 


2 


1 


2 


1 


2 


4 




Pa 


1 


2 


1 


2 


1 


2 


1 


2 


4 


2 




Pi 


1 






1 




1 


1 




1 


1 




P 2 


2 


1 


1 


2 


1 


2 


2 


1 


3 


3 




Pa 


1 


2 


2 


1 


2 


1 


1 


2 


3 


3 



uniquely determine the symmetry of the structure, inde- 
pendently of the kind of the property taken into account. 
The type of the phase transition and the property un- 
der consideration are encoded in the form of the basis 
vectors. 

In this work the symmetry analysis is applied to the 
structure of the cr-phase described by the P42/mnm space 
group, where the Wyckoff positions 2a, 4/, 8i, 8V and 8j 
are occupied by both Cr and Fe atoms. The aim of our 
analysis is to find a possible ordering of the magnetic mo- 
ments that allow for a magnetic phase transition without 
any change of the crystal structure. As the first step, the 
active representations of the invariable lattice k= (0,0,0) 
for different positions and for different types of ordering 
(modes) have been calculated and shown in Table IIVI 
The designations Pi, P2 and P3 of the mode types corre- 
spond to a scalar (change of the probability of the site oc- 
cupations), polar (displacement of atoms from the equi- 
librium positions in the high symmetry structure) and 
axial (ordering of the magnetic moments) modes, respec- 
tively. Hence, for a discussion of a possible magnetic 
structure in the cr-phase, only the latter modes are rele- 
vant. The values presented in Table lTVl indicate a number 
of different basis vector sets available for a given IR. Since 
the ordering of the magnetic moments is not associated 
with any change of the crystal structure in this system, 
only the models of the magnetic structure belonging to 
T\ IR are considered. 

As can be seen from Table HVl the ordering of the mag- 
netic moments in this case is possible only on positions 
8i (C and D sites) and 8j (E sites). The results of our 
calculations are given below: 

Site 8i: 

Hi = (0, 0, fi ) = fi 2 = M3 = ^4 = -Ms = -Me = ~M7 = 
-Ms 
Site 8j: 

Mi = (Mo, -Mo, 0) = -M2 = ~M7 = Ms 
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FIG. 5: (Online color) Magnetic moments of Fe and Cr atoms for five crystallographic sites versus the number of NNFe atoms. 
Circles and solid lines stand for the KKR results (this work), whereas triangles and diamonds correspond to the calculations 
reported in Refs. andQj], respectively. 
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FIG. 6: (Online color) Magnetic moments of Fe and Cr atoms for five crystallographic sites versus the number of NNf e atoms 
as obtained with the APM model of the magnetic ordering (circles). Solid lines stand for the best fits to the magnetic moments 
in the FM ordering shown for the sake of comparison. 



= (Mo, Mo, 0) = -m = n 6 = -// 5 

where the numbering of the atomic positions is 
given according to the International Tables for X-Ray 
Crystallography^. The parameter /zo is free (both in 
magnitude and in sign) what means that the symme- 
try analysis is not able to predict differences neither in 
values nor in relative ordering of the magnetic moments 
between different sublattices. Due to this freedom differ- 
ent models of possible magnetic orderings may be con- 
sidered. However, the calculations clearly show that the 
arrangements of the magnetic moments within 8i and 



8j sublattices are antiferromagnetic with strictly given 
sequences, and, additionally, at 8j site the structure is 
noncollinear (4 magnetic moments aligned antiparallel in 
the basic plane along the diagonal and 4 other ones also 
antiparallel but perpendicular to the previous moments) . 
It is worth noticing that the purely antiferromagnetic 
state suggested by the aforementioned symmetry analysis 
is possible when atomic positions are occupied by given 
types of atoms with the same probability. This type of 
the lattice occupation is broken in the case of a Fe/Cr dis- 
ordered alloy, and the proposed magnetic structures can 



be treated as initial conditions for the electronic structure 
calculations. This is the reason why, one of the consid- 
ered model of the magnetic structure, has been called 
APM, instead of antiferromagnetic one. 



D. Ordered model of magnetic APM a-FeCr 

The KKR calculations using the APM model were per- 
formed assuming the input crystal potentials obtained 
for the FM state. The computations started with ini- 
tial spin-polarized potentials, but this fact had no effect 
on the finally converged APM results. Such procedure 
allowed us to substantially accelerate the self-consistent 
convergence. According to the APM model, antiferro- 
magnetic ordering was assumed on the sublattices C and 
D. The symmetry analysis does not provide any infor- 
mation about the relative magnetic moments sequences. 
These sequences can be known within each sublattice sep- 
arately. Consequently, two separate magnetic structures 
should be considered with the following relative sublat- 
tice settings: C (till) and D (tltl) as well as C (tltl) 
and D (4-tit) > which led us to different values of the mag- 
netic moments. Since only collinear arrangement was al- 
lowed in the KKR calculations, antiparallel coupling of 
the magnetic moments on E site was not taken into ac- 
count. 

As is clear from Fig. [6l the assumed APM model for 
the Fe/Cr magnetic moment arrangement has converged 
to a stable magnetic configuration with atomic moments 
of the opposite polarization, which can be especially well 
seen in the case of Fe atoms located on sublattices C and 
D (and also to some extent on sublattice E). Furthermore, 
the APM model computations have generated interme- 
diate Fe magnetic moments with a different alignment, 
particularly on sublattices A and E, which were not ob- 
served for the FM ordering. There was not revealed any 
visible change on the sublattice B with respect to FM. 
The changes of the magnetic moments of Cr atoms on all 
sublattices have been found similar to the Fe ones. The 
smallest effect has been observed on the sublattice E. 

In the case of APM model, the total magnetic mo- 
ment as determined by averaging the moments of indi- 
vidual atoms yielded 0. 21 fiB, with negligible contribu- 
tion (c.a. — 0.01/xb) from the outside muffin-tin volume, 
which ultimately gives 0.20/is (Table |TTT}. This value 
is significantly smaller than that obtained with the FM 
model (0.59/is), and much closer to the experimental 
value (0.14/xs)^. This suggests that the APM-ordering 
may be responsible for the small value of the magnetic 
moment deduced from the magnetization measurements. 

A spatial distribution of the magnetic moments in the 
a— phase unit cell found with the APM model shows a 
behavior similar to that detected with the FM model 
(see Fig. 2b). In this case, however, the projected areas 
of the magnetic moments with the same orientation are 
markedly smaller. 
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FIG. 7: (Online color) Total electronic DOS in the 
a— Fei6Cri4 alloy obtained with the FM (blue, dashed) and 
APM (black, solid) models, as calculated by the KKR-CPA 
method. 



E. Chemically disordered model for the APM 
ordering 

Bearing in mind the aforementioned KKR results both 
in the FM and the APM ground states, we have also stud- 
ied the effect of a chemical disorder on electronic struc- 
ture and magnetic properties of the cr— Fei6Cri4 alloy 
using the KKR-CPA method. 

In this case, the disorder (treated as random) was 
introduced on three sublattices (B: Feo.25Cro.75, C: 
Feo.25Cro.75 and E: Feo.375Cro.625), where the disorder is 
the most significant as follows from the neutron diffrac- 
tion experiment (see Table [D} . Due to the highly time- 
consuming calculations, the sublattices A and D were left 
fully occupied by Fe atoms. This modification required 
to slightly decrease the Fe concentration on other sublat- 
tices, in order to maintain the FeigCri4 stoichiometry for 
the direct comparison between the KKR and the KKR- 
CPA results. The KKR-CPA results are discussed here 
for the FM and the APM models, neglecting the non-spin 
polarized computations. 

The total DOSs computed for the cr-Fei 6 Cri 4 (Fig. [7]) 
are quite similar in both cases. However, there are few 
subtle differences between the two presented DOS spec- 
tra. In the FM state, the DOS is more spin polarized (up- 
DOS and down-DOS arc slightly shifted into lower and 
higher energy range, respectively) as well as two down- 
DOS peaks are more pronounced. This feature yields the 
total magnetic moment as large as ~ 0.4/is (~ 12.0/is 
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per unit cell) against only ~ 0.15/is (~ AA^b per unit 
cell) in the APM state (TableHn]). Interestingly, the DOS 
at Ei? in the APM state is a bit smaller than that in the 
FM one, which may tentatively support a better stability 
of the APM state, relative to the FM one. Before inspect- 
ing in more detail the local origin of the magnetism in the 
disordered a— FeigCri4 alloy, let us remark that disorder 
effects markedly decrease the magnetic moments on Fe 
atoms on all sublattices, and, in consequence, the total 
magnetization per unit cell. This feature is well seen for 
the FM state, where the value of OA/is was calculated 
with KKR-CPA for the total magnetization, compared 
to almost 0.6/is with KKR as obtained via averaging 
over a set of results obtained for the ordered approxi- 
mant approach. Quite similar behavior was found for the 
APM state, since the magnetization values gained from 
the KKR-CPA and the KKR methods also differ (0.15/is 
for disordered alloy against 0.20 [1b for ordered models 
of the cr— FeigCri4). Hence, an inclusion of the chemical 
disorder into the electronic structure calculations has ad- 
ditionally lowered the total magnetization by about 25%. 

The choice of the APM magnetic structure for the 
KKR-CPA calculations has resulted in the total magnetic 
moment (0.15 hb) in an excellent agreement with the ex- 
perimental data (0.14 /xb) for the a— Feo.53sCro.462 alloy. 
It should be added that in the APM model, according to 
the symmetry analysis, sites C and D with multiplicity 8 
were split into two "sub-sites", each with multiplicity 4, 
having an antiparallel orientation. Consequently, in or- 
der to distinguish the magnetic moments and DOSs of Fe 
atoms on D site and those of Fe/Cr atoms on C site that 
are arranged parallel or antiparallel to other sites, (A, B 
and E), we have introduced indexes u and d, respectively. 

Inspecting the Fe DOSs on all sites in the APM state 
(Fig. |5]), one can notice two sites (B and C) with strongly 
polarized electronic spectrum, giving rise to large local 
magnetic moments of +2.08 /is (B site) as well as of - 
2.06 hb (Fed) and of +2.11 /z B ( Fe„) on C site. The 
spin polarization of DOS on other sites remains much 
smaller, resulting in relatively lower values of Fe mag- 
netic moments on A (+0.32 /is), D (-0.44 on Fe^ and 
+0.65 on Fe M ) and E (+0.90 fi B ). The KKR-CPA 
results well reflect the aforementioned correlation (Sec. 
Ill B) between the values of the Fe magnetic moment 
and the average distance to NN atoms, well supporting 
the KKR findings for the ordered approximants. It is 
worth noticing that on these sites, where Cr atoms were 
substituted (B, C and E), a small Cr magnetic moment 
was found to be always aligned antiparallel to the Fe mag- 
netic moment on the same site. Such electronic structure 
behavior was also found for the FM state using the KKR- 
CPA method. As already mentioned, the polarization of 
the Cr DOS is weak and for E site a slight DOS asymme- 
try for two spin directions yielded the magnetic moment 
of about -0.20 fiB- In the case of B and C sites, the cal- 
culated Cr magnetic moment is as small as -0.10 /is, on 
average. 

At present, we must rely on the above described inter- 
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FIG. 8: (Online color) Site-decomposed KKR-CPA DOS for 
the a— Fei6Cri4 alloy as obtained with the APM model. 



pretation of the magnetic ordering in the a— FeCr sys- 
tem, since first the APM structure appears the simplest 
one and second there are no available experimental evi- 
dences indicating more complex (e.g. noncollinear) mag- 
netic structure. 



IV. CONCLUSIONS 

The results of the electronic structure calculations us- 
ing KKR and KKR-CPA methods reported in the present 
study can be summarized. Both Fe and Cr atoms in 
the a— Fei6Cri4 alloy are magnetic. Magnetic moments 
of Fe atoms for each lattice site linearly increase with 
the number of Fe atoms in the nearest-neighbor shell, 
NNi? e , whereas Cr atoms are coupled antiparallel with 
respect to Fe atoms (values of the magnetic moments for 
both Fe and Cr atoms are characteristic of a given lat- 



10 



tice site). Ferromagnetic ordering of the Fe-site magnetic 
moments accounted for the calculations results in a sig- 
nificant overestimation of the magnetism in the investi- 
gated sample. Fortunately, admission of an anti-parallel 
ordering on some selected sites, suggested by the sym- 
metry analysis, in a combination with a partial chemical 
disorder yields the excellent agreement with the exper- 
imentally found value of the average magnetic moment 
per atom. 
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